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ABSTRACT 

We present results from a high-resolution and large-scale hybrid (fluid electrons and particle-in-cell protons) 
two-dimensional numerical simulation of decaying turbulence. Two distinct spectral regions (separated by a 
smooth break at proton scales) develop with clear power-law scaling, each one occupying about a decade in 
wave numbers. The simulation results exhibit simultaneously several properties of the observed solar wind fluc¬ 
tuations: spectral indices of the magnetic, kinetic, and residual energy spectra in the magneto-hydrodynamic 
(MHD) inertial range along with a flattening of the electric field spectrum, an increase in magnetic com¬ 
pressibility, and a strong coupling of the cascade with the density and the parallel component of the magnetic 
fluctuations at sub-proton scales. Our findings support the interpretation that in the solar wind large-scale MHD 
fluctuations naturally evolve beyond proton scales into a turbulent regime that is governed by the generalized 
Ohm’s law. 

Subject headings: The Sun, Solar wind, Magneto-hydrodynamics (MHD), Plasma, Turbulence. 


1. INTRODUCTION 

In-situ measurements of the solar wind plasma and electro¬ 
magnetic field show spectra with a power-law scaling_span- 
ning s everal decades in frequency, f (e.g. lAlexandrova et al.1 
2009: ISahraoui et al~]l201Ct (Roberta 1201 Oh . Power-laws sup¬ 
port an interpretation in term of turbulent fluctuations, al¬ 
though the rich variety of spectral features is not easily ex¬ 
plained in the framework of known turbulent theories and phe- 
nomenologies. 

For frequencies in the so-called magneto-hydrodynamic 
(MHD) range, 10 -4 Hz < / < 10 -2 Hz at 1 AU, the magnetic 
field spectrum and the kinetic field spectrum have a different 
scaling, t he former being proportional to f ~ 5/3 while the latter 
to f~ 3/2 dPodestae t al. 2007]; ISalem et all 120091: IWicks et aD 
1201 lk iTbssein et al.ll2009l) . While a magnetic excess is gener¬ 
ally found in solar wind turbulence, only recently the spec¬ 
trum of residual energy (the difference between magnetic 
and kinetic energy) was shown to have a power-law scaling 
with a spectral index -2 (IChen et aOl2013al) . Such finding 
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confirms early predictions on the residual energy spectrum 
( Gra ppin et al JI1983J) and numerical results o f incompressible 
MHD simulations (IMuller & Grappiiil 120051) . Note that the 
three spectral indices (-3/2, - 5/3, - 2) for the kinetic, 
magnetic, and residual energy spectrum are not reproduced si¬ 
multaneously in any direct numerical simu lation (DNS) (e.g. 
IMuller & Grappinll2005l: IChen et al.ll201 lbh unless a particu - 
lar driving is applied to large scales dBoldvrev et~aLll2011[) . 
Finally, in the MHD range, magnetic and velocity fluctua¬ 
tions are dominated by the transverse components with re¬ 
spect to the ambie nt magnetic field B o (e.g. Smith et al. 2006; 
IWicks et alil20Tl . 

Moving to higher frequencies, / > 10 -2 Hz, there is grow¬ 
ing evidence that kinetic effects become important and change 
the nature of the self-similar spectra of fluctuations observed 
for / > 1 Hz. A spectral break appears in magnetic and 
velocity spectra at proton scales, separating the MHD iner¬ 
tial range cascade from a second power-law interval at ki¬ 
netic scales. The physical scale associated with this spectral 
break has not been iden t ified vet (e.g.lBourouaine et al.ll2012l: 
iBruno & Trenchill2014l: IChen et al .1120141) . The spectral in¬ 
dex of m agnetic fluctuations after the break v aries between 
(-4, -2) dLeamon et al.lll998l: ISmith et aPl2006l) . although it 
tends to cluster around a slop e of -2.8 for higher frequen¬ 
cies lAlexandr ova et al.ll2012l) . The change in the t urbulence 
regimes also shows up in the density spectrum ( Che n et al.1 
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2013blh which steepens and couples to the parallel compo¬ 
nent of magnetic field. The latter becomes as energetic as the 
two perpendicular components, resul ting in an increase of the 
so-called m agneti c compressibilitv jAlexandrova et al.ll2008i: 


ISalem et aPl2012l: iKivani et aDl2Q13h . Finally, measurements 

at 1 AU show that the spect r um of electric field flattens at 
about 1 Hz dBale et al.ll2QQ5b [Kellogg et aD 120061) . although 
the noise level hinders the determination of a precise spectral 
scaling. 

The m easure of structure funct i ons of third order at MHD 
scales (ISorriso-Valvo et al.l [2 007: MacBri de et all l200~8h and 
of hig h-order at MH D (Sal em et all120091) and at sub-proton 
scales (IKivani et al.ll2Q13h provided additional evidences that 
fluctuations are turbulent all the way down to electron scales 
in the solar wind. While DNS are able to reproduce some 
aspec ts of either the MHD r ange (e.g. [M aron & Gold reichl 

2001; Mason et al.l 120081: iBeresnvak & Lazarian 


Grappin & Muller 2010: Leeetal . 2010; Chen et al. 12011 


Bold yrev et all 1201 HT Dong et all 12014 ) or the sub - proton 


Markovskii & Vasauez 

2011 : 

Camporeale & Burgess 

2011 ; 

Boldvrev et all 2013 

J 

GarvetalJ 2012 

; Wan et al. 

20121 


Ml, to our knowledge a clear indication that a turbulent 


regime establishes in the whole spectrum spanning the two 
ranges has never been reported so far. 

In this work we present results from a high-resolution 
hybrid (fluid electrons, particle-in-cell protons) two- 
dimensional (2D) DNS of turbulence and provide the first 
direct numerical evidence of the simultaneous occurrence of 
several features observed in the solar wind spectra. These 
include i) the different scaling of magnetic and kinetic 
fluctuations in the MHD range, ii) a magnetic spectrum with 
a clear double power-law scaling separated by a break, iii) 
an increase in magnetic compressibility at small scales, iv) a 
strong coupling between density and magnetic fluctuations 
at small scales. The electric field spectrum is also consistent 
with observations, showing a change in the spectral properties 
at sub-proton scales. Our results indicate that the switch in 
the spectral slopes observed in the solar wind results from 
the natural continuation of a large-scale MHD turbulent 
cascade through proton and down to electron scales, where 
the different field couplings are governed by the non-ideal 
terms of the Ohm’s law. 


2 . NUMERICAL SETUP 

The kinetic model uses the hybrid approximation: electrons 
are considered as a massless, charge neutralizing, isother¬ 
mal fluid; ions are described by a particle-in-cell model (see 
[Matthews! 19941 for detailed model equations). The character¬ 
istic spatial and temporal units used in the model are the pro¬ 
ton inertial length d p = va/Q p , va being the Alfven speed, and 
the inverse proton gyrofrequency 1/Q p , respectively. We use 
a spatial resolution Av = Ay = 0.125 d p , and there are 8000 
particles-per-cell (ppc) representing protons. The resistive co¬ 
efficient is set to the value 77 = 5 10 -4 4 nvAC~ l fl~ l to prevent 
the accumulation of magnetic energy at the smallest scales. 
Fields and moments are defined on a 2D x-y grid with dimen¬ 
sions 2048 2 with periodic boundary conditions. Protons are 
advanced with a time step At = 0.025 Q" 1 , while the magnetic 
field B is advanced with a smaller time step At B = At/ 10. 
The number density n is assumed to be equal for protons and 
electrons, n p = n e = n, and both protons and electrons are 



Fig. 1.— Root-mean-square parallel current density (black line), perpen¬ 
dicular and parallel magnetic fluctuations (red-solid line and red-dashed lines 
respectively), and density fluctuations (purple line) as a function of time (nor¬ 
malized to the inverse of the proton gyrofrequency Q p ). As a reference, the 
nonlinear time at the initial time is about 20 Q~ l . 


isotropic, with J3 p = J3 e = 0.5 where p p , e = %7inK B T p , e /are 
the proton (electron) betas (here K B is the Boltzmann’s con¬ 
stant, #0 the ambient magnetic field, and T p , e are the proton 
and electron temperatures). 

We impose an initial ambient magnetic field Z?o = Bqz, per¬ 
pendicular to the simulation plane. We add an initial spec¬ 
trum of linearly polarized magnetic and bulk velocity fluctu¬ 
ations u, with only in-plane components. Fourier modes of 
equal amplitude and random phases are excited in the range 
-0.2 < k Xt y < 0.2, assuring energy equipartition and vanishing 
correlation between kinetic and magnetic fluctuations. Ini¬ 
tial velocity fluctuations have vanishing divergence and den¬ 
sity fluctuations are also vanishing (in the limit of numerical 
noise). Quantities are defined as parallel (||) and perpendic¬ 
ular (_L) with respect to Bq. We define the omnidirectional 
spectra, 

EAk±) = 6'¥\k ± )/k ± = E V 2 2 D (k±l (1) 

I*lI =*l 

where ¥ are the Fourier coefficients of a given quantity 'F 
(we use E and J to indicate electric field and current density 
respectively) and dT^A^) is the amplitude of the fluctuation 
¥ at the scale k_ L . We also define the root mean square value 
(rms) as 

x ¥ rms = V0? 2 > - ('P> 2 - (2) 

where (...) stands for a real-space average over the whole sim¬ 
ulation domain. With these definitions, the initial conditions 
have E u ~ E B cc k_ L with B rms /Bo ~ 0.24 allowing for a fast 
turbulent dynamics sustained for about 300 D" 1 (the nominal 
nonlinear time at the beginning of the simulation is approxi¬ 
mately 20 Q” 1 , but it increases at later time). 

3. RESULTS 

In Figure |T] we plot the rms of the parallel current den¬ 
sity, of the parallel and perpendicular magnetic field, and of 
the density fluctuations. The current density increases until 
t = 200 Q" 1 , reflecting the formation of small scales due 
to the development of a turbulent cascade, and then declines 
smoothly. The decay is slow, since larger and larger scales 
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Fig. 2.— Contour plot of the perpendicular magnetic energy at t - 200 . 


continue to feed the cascade at later times. Accordingly the 
perpendicular magnetic field declines steadily after a transient 
increase. Shortly after the beginning, fluctuations in the par¬ 
allel component of magnetic field and in the density appear, 
slowly increase, reaching a shallow maximum at the same 
time of the current density, and then decline slowly. The initial 
growth is due to the generation of a low level of compressive 
fluctuations. Velocity fluctuations (not shown) behave simi¬ 
larly to magnetic fluctuations, with the perpendicular compo¬ 
nent declining monotonically (there is no initial growth) and 
the parallel component originating from compressive effects. 
In the following we will show spectra at the time of the peak 
of the current density t - 200 Q" 1 , but all the turbulent proper¬ 
ties are stable and remain valid until the end of the simulation 
(t = 500 Q; 1 ). 

In Figure [2] we show isocontours of the perpendicular mag¬ 
netic field energy. This snapshot highlights the formation of 
intense vortex-like and filamentary structures. The latter re¬ 
flect the local anisotropy of small scales fluctuations, while 
their random orientation assures the statistical isotropy of the 
two-dimensional spectrum: we thus consider in the following 
only omnidirectional spectra. 

In Figure [3] (top panel), we show the spectra of the total ve¬ 
locity, magnetic, and electric field. The magnetic spectrum 
(red line) has a double power-law scaling, each power-law 
range occupying about one decade, with a break at k±d p ~ 2 
that separates the MHD from the sub-proton range. The bulk 
velocity spectrum (blue line) also has a power-law scaling in 
the MHD range but it falls off abruptly at k ± d p ~ 1, not show¬ 
ing any clear power-law at higher wavenumbers. At smaller 
scales it reaches the ppc noise level, estimated as the level of 
velocity fluctuations at t = 0 (light blue dashed line). Finally 
the electric field spectrum (green line) follows the velocity in 
the MHD range (k ± d p < 0.4) and tends to flatten as it enters 
the sub-proton range (k ± d p > 2). 

These spectral properties are qualitatively and quantita¬ 
tively in agreement with observed solar wind spectra. In the 
MHD range the magnetic and kinetic spectra are power-laws 
with scaling consistent with Eb oc k~ ± 5/3 and E u a: k ± 3/2 , re- 
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k±d p 


Fig. 3. — Top panel. Omnidirectional spectra of total magnetic (red), total 
kinetic (blue), and total electric field (green) fluctuations versus perpendicu¬ 
lar wavenumber k± at t = 200 Q" 1 . The spectrum of kinetic energy at t = 0 is 
plotted with a dashed light-blue line as an indicator of ppc noise level.Dashed 
black lines are references for the corresponding spectral indices. Bottom 
panel. Magnetic (red), kinetic (blue) and residual energy (black) spectra com¬ 
pensated by k 5 / 3 , k 3 / 2 and k]_, respectively. 


spectively, as can be seen in the bottom panel of Figure [3] 
where the spectra are compensated by k 5/3 and k 3/2 respec¬ 
tively. In the same panel we also plot the residual energy spec¬ 
trum, Er = Eb~E u , which has a power-law scaling over about 
one decade in the MHD range with a spectral index « -2 as 
in observations dChen et al.lll2013al) . In addition, in the sub¬ 
proton range the magnetic spectrum scales as E B oc k “ 3 , a 
spectral index which is very close to the value -2.8 reported in 
observations ( Alexandr ova et alJl2009h . Note that the electric 
field spectrum is strongly coupled to the bulk velocity spec¬ 
trum at MHD scales (they are basically indistinguishable for 
k±d p < 0.4), reflecting the dominance of the ideal MHD term 
(I u X B\ ~ B 0 u ± ) in the generalize d Ohm’s law, and con¬ 
sistent with solar wind observations dChen et all 1201 la ). At 
smaller scales, it decouples from the velocity spectrum since 
the Hall term (/ x B/n) and the electron pressure gradient 
term (VP e /n) start to dominate. 

Since both other fields and derivatives enter in its compu¬ 
tation, E is the field that is mostly affected by numerical ef¬ 
fects and it’s not straightforward to give a simple estimate of 
its noise level, as done for the velocity field. Ultimately, we 
can reasonably claim that the shallower slope of its spectrum 
for 2 < kj_d p < 7 is of physical nature, while its behavior at 
smaller scales is most likely not. On the contrary, quantitative 
results for the spectra of magnetic and density fluctuations are 
more robust even at larger wave numbers. A detailed descrip¬ 
tion and discussion about different sources of numerical noise, 
e.g. the finite number of ppc, will be given in a companion pa- 
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Fig. 4.— Omnidirectional spectra of density (purple line), parallel magnetic 
field (red line), parallel and perpendicular electric field (green dashed and 
solid line, respectively) versus perpendicular wavenumber k± at t = 200 Q~ l . 


per dFranci et al .1120 151) . For the purpose of this Letter, what 
matters is that such numerical noise does not affect either the 
qualitative scaling of the electric field spectrum for k ± d p < 1 
or the estimate of the spectral indices of other fields up to 
k ± d p ~ 10 (except the velocity field which is presumably af¬ 
fected by the noise level at k ± d p > 4). 

The transition from the MHD regime to the sub-proton 
regime is not only characterized by a change in the spectral 
indices, but also by an increase of energy of the parallel mag¬ 
netic field and the density fluctuations relative to other fields. 
These are shown in Figure 01 along with the parallel and the 
perpendicular electric field spectrum. The density and paral¬ 
lel magnetic fluctuations are coupled in the whole range of 
scales. In the MHD range, they have a flat spectrum that is 
an order of magnitude smaller than the perpendicular elec¬ 
tric field. This also results in a small power in the spectrum 
of the total magnetic field intensity E\ B \ < E B (not shown), 
consi stently with solar wind observations ( Horbury & Balogh 
1200ll) . In the sub-proton range, E B \\ and E n steepen, both hav¬ 
ing a clear power-law scaling with index -2.8. By comparing 
Figures [ 3 ] 0 ] one can see that the parallel and perpendicular 
components of magnetic fluctuations become comparable at 
the sub-proton scales, leading also to E\ B \ ~ E B . Concern¬ 
ing the electric field spectrum, at all scales the perpendicu¬ 
lar component E e± dominates by a factor ~ 100 the parallel 
component E E \\, reflecting the fact that in our configuration 
the leading terms of the generalized Ohm’s law are linear and 
quadratic in the fluctuations’ amplitude for E e± and E E \\ re¬ 
spectively. Note that E e ± flattens at the sub-proton scales 
and Ep\\ steepe ns in qualitative agreement with observations 
(iMozer & ChenllMl3h . It is hard to determine the spectral 
index of E e± at sub-proton scales; a rough estimate gives 
oc k~ 0 - 8 , consistent with E being determined by the Hall and 
pressure terms. In fact, retaining only the leading order in the 
expression of E one gets E e ~ E e± oc k\E B ^ >n ~ k“ 0,8 . 

We can further compare our results with observations con¬ 
sidering three non-dimensional ratios involving density, mag¬ 
netic and electric field fluctuations shown in Figure 0 Con¬ 
sider first the magnetic compressibility, the ratio of parallel 
to total magnetic fluctuations (red line). It is negligible in 
the MHD range, increases while approaching the sub-proton 
scales, and finally saturates to a level 6B\\/6B ~ 0.5. Thus, 
magnetic fluctuations have mainly perpendicular components 



kj_d p 


Fig. 5. — Ratio of perpendicular electric field to perpendicular magnetic field 
(green), ratio between normalized density and perpendicular magnetic fluc¬ 
tuations (purple, see text for the normalizations), and ratio of parallel to total 
magnetic field fluctuations (magnetic compressibility, red). The numerical 
noise affects the ratios for kj_d p >l (vertical dotted line). 


in the MHD range but tend to become isotropic at small 
scales, approaching a value SBj ~ 6B]_/3 , which is within 
the range (~ 0.2-^0.5) measured in the solar wind at space¬ 
craft frequencies larger then 1 Hz dKivani et al.ll2013l) . This is 
also in very good agreement with the level of magnetic com¬ 
pressibility expected for kinetic Alfven wave turbulence for 
the parameters adopted in our simulation (e.g. lBoldvrev et al.l 
2013). 

The purple line in Figure \5\ shows the ratio of nor¬ 
malized density fluctuations over normalized perpendicu¬ 
lar magnetic fluctuations, Sh/SB ± , where SB ± = SB ± /Bo 
and 6h = TSn/n 0 respectively, and T (3/4 in our simula¬ 
tion) is a non-dime nsional kinetic normalizat i on that depends 
on T p,T e ,/3 p ,VA (ISchekochihin et al.1120091: iBoldvrev et all 
120131) . With this normalization SB ± and Sn are expected to 
have the same amplitude for kinetic Alfvenic fluctuations. In¬ 
deed Sh/SB ± increases and then saturates at a value ~ 1 at 
sub-proton scales. Note that the plateau and its value ~ 1 are 
consistent with observations (on average Sh/SB ± = 0.75, cf. 
Chen et all 2 013b ). 

Finally, we plot the ratio between the perpendicular electric 
fluctuations (normalized by the Alfven speed) and the perpen¬ 
dicular magnetic fluctuations (green line). Similarly to the ob¬ 
served frequency spectra in the solar wind frame dBale et al.l 
2 0051) . this ratio is about 1 in the MHD range, where the MHD 
term (u x B) dominate. At k ± d p ~ 1 the ratio increases re¬ 
flecting the role of the Hall term (J x B/n) and the pressure 
gradient term (VP e /n) in the generalized Ohm’s law. 

4. CONCLUSION 

In this Letter we show that hybrid 2D large-scale, high- 
resolution simulations of turbulence are able to reproduce si¬ 
multaneously several aspects of the MHD range and of the 
sub-proton range of solar wind spectra. 

Two noticeable examples are given by the spectra of the 
magnetic field and of the electric field. The former displays 
a clear double power-law scaling, with spectral indices —5/3 
and -3 in the MHD and sub-proton range respectively, sepa¬ 
rated by a smooth break at k ± d p ~ 2. The electric field spec¬ 
trum also shows a change in the spectral properties at about 
the same scales, being coupled to velocity fluctuations in the 
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MHD range, and becomes shallower at sub-proton scale. It 
is also worth noting that in the MHD range we found the 
scaling observed in the solar wind for the magnetic, kinetic, 
and residual energy spectra (respectively -5/3, - 3/2, and 
-2). To our knowledge this is the first time that these spec¬ 
tral indices are obtained for turbulence with vanishing cor¬ 
relation between magnetic and velocity fields. DNS of in¬ 
compressible MHD usually capture only the scaling of the 
residual energy and the total energy dMuller & Grappinl2005l) 
while Reduced MHD fails in reproducing velocity and kinetic 
spectral indices (IChen et al.ll201 lbh or requires special driving 
dBoldvrev et al.11201 ll) . This may indicate that it is necessary 
to go beyond the incompressible MHD approximation even in 
the inertial range. Further work is needed to test this possibil¬ 
ity, extending the analysis to a full 3D simulation. 

In the sub-proton scales we found an increase in magnetic 
compressibility and a strong coupling between density and the 
parallel component of magnetic fluctuations - both having the 
same spectral index of -2.8 - with the main cascade of E B± 
driven from the MHD scales. All these spectral indices match 
or are consistent with observations. The only relevant dis¬ 
crepancies are the flat spectra (slope ~ 0) of parallel magnetic 
fluctuations and density fluctuations in the MHD range. In the 
solar wind they have a spectral index -5/3 (e.g. IChen et al.l 
2 012\ ). This aspect is not fully captured by our simulations 
probably because of the limited compressibility imposed by 
the 2D dynamics and/or by the value of the proton Note 
however that this does not prevent the full development of 
a compressible cascade at kinetic scales, in good agreement 
with observations. 

Properties shown in Figure [5] are consistent with the tur¬ 
bulence at sub-proton scales being ruled by fluctuations with 
properties of kinetic Alfven waves. However, note that the 
level of magnetic and gas compressibility expected for this 


regime follows from more general properties of the thermody¬ 
namical state assumed for the plasma (J3 e , and ion-electron 
temperature ratio), which govern the couplings between the 
different fields B , E and n via the generalized Ohm’s law. 
In the low-frequency regime (i.e. below the whistler range), 
the ratios ; Sh 2 /SB ' 2 ., SB\\/6B are not expected to depend on 
k dBoldvrev et al.l 120131) since they do not rely on the spe¬ 
cific dispersion relation of the fluctuations. In this sense, the 
plateaus at 2 < k± < 7 in Figure [5] represent a more gen¬ 
eral and likely universal manifestation of low-frequency tur¬ 
bulence at kinetic scale, and this is how we intend to present 
them here. 

As a concluding remark, we stress that our simulation im¬ 
plements a finite resistivity to assure a source of damping at 
small scales for the magnetic fluctuations, and thus to prevent 
energy accumulation and the consequent artificial flattening 
of the spectrum. Although a more detailed and quantitative 
analysis of the related effects will be given in a forthcoming 
paper dFranci et^0l2Q15h . we anticipate that the values of re¬ 
sistivity and the number of ppc affect the ion heating proper¬ 
ties. 
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